UUITP-03/13 



Fully stable dS vacua from generalised fluxes 

Johan Blaback, Ulf Danielsson and Giuseppe Dibitetto 

Institutionen for fysik och astronomi, 

Uppsala Universitet, 
Box 803, SE-751 08 Uppsala, Sweden 
{johan. blaback, ulf . danielsson, giuseppe .dibitetto}@physics .uu. se 



ABSTRACT 

We investigate the possible existence of (meta-)stable de Sitter vacua within 
M = 1 compactifications with generalised fluxes. With the aid of an algorithm 
inspired by the method of differential evolution, we were able to find three 
novel examples of completely tachyon-free de Sitter extrema in a non-isotropic 
type IIB model with non-geometric fluxes. We also analyse the surroundings 
of the aforementioned points in parameter space and chart the corresponding 
stability regions. These happen to occur at small values of the cosmological 
constant compared to the AdS scale. 



Introduction 



The possible existence of (met a-) stable de Sitter (dS) vacua in string theory, is considered 
to be of crucial importance if we want to achieve agreement with the main cosmological 
observations of the last two decades. The problem has been addressed following various 
approaches using string compactifications, and their link with the corresponding effective 
descriptions. Focusing on particular effective supergravity descriptions, which often are 
referred to as STU-models, the introduction of fluxes turns out to be the key, at least 
perturbatively, for inducing a superpotential of the right form [TJ. 

Within these supergravity models, several theories have been studied that describe string 
backgrounds in type IIA (2 13 , and type IIB 14 19 . So far, there are no known example 
of a fully meta stable dS vacua. All geometric compactifications have unstable directions 
already within the isotropic truncation of the STU-model. Using non-geometric fluxes (20l 
it was possible to construct the first example of an isotropically stable dS extremum [18] . 
This solution is still, however, unstable in the non-isotropic directions. This situation was 
recently analysed in a more systematic way in ref. [2~T] . 

The main result of this paper are the first examples of fully stable dS vacua in the STU- 
model. We show that non-geometric fluxes provide enough freedom in parameter space 
for tuning the masses of all the moduli to be positive, even the non-isotropic ones. In 
appendix [A] we analyse the whole region of stable dS in parameter space around such stable 
dS points by providing some insightful plots. In appendix [Bj we give some details concerning 
the algorithm that allowed us to find the aforementioned stable dS vacua. 



The 7j2 x Z2 orbifold with generalised fluxes 

Orbifold compactifications of type IIB string theory on T 6 /(Z2 x Z2) with 03/07-planes 
(and duals thereof) with generalised fluxes can all be placed within the same framework of 
effective four-dimensional supergravity descriptions that are known as STU-models. These 
theories enjoy M = 1 supersymmetry and SL(2) 7 global bosonic symmetry. The action of 
such a global symmetry on the fields and couplings can be interpreted as the effect of string 
dualities. 

The scalar sector contains seven complex fields spanning the coset space (SL(2)/SO(2)) 7 
which we denote by $ Q = (S,Ti,Ui) with i = 1,2,3. The kinetic Lagrangian follows from 
the Kahler potential 

3 3 
K = - log (-i (S - S)) - J2 lo S H ( T * - T $ ~ E lo S H & - Ui)) . (1) 



This yields 

{-i{S -S)) 2 U\ (-i(T, - Tij) (-i(U t - I/,)) 2 / 

The presence of fluxes induces a scalar potential V for the moduli fields which is given 
in terms of the above Kahler potential and a holomorphic superpotential W by 



V = e K (-3|WT + K ap D a WDpW) , (3) 

where K af5 is the inverse Kahler metric and D a denotes the Kahler-covariant derivative. 

The general form of a superpotential induced by locally geometric fluxes in type IIB with 
03/07-planes is given by 

W = P F (Ui) +S PuiUi) + J2 T k Pq\ud + S Y,T k P { p ) {U l ) , (4) 

F flux H flux h Q fl ux k P flux 

where Pp, Ph, Pq^ and P p k ^ are cubic polynomials in the complex structure moduli given 
by 

p F {Ui) = a -Y,^ ) u i + Y,4 )J ^^ 1 - ^u,u 2 u, 
p H {Ui) = -b, + Y.bfu l -Y.bf l ^^ + hu 1 u 2 u, 

m = of + e u< - e cf ) - 4 fc) «/i ^ c/3 , (5) 

i i U% 

ppm = -df - e df > ^ + e rff } + 4 fe) t/i ^ ^ • 

The IIA and IIB interpretation of the above superpotential couplings is summarised in 
table |U 



Fully stable de Sitter vacua 

The setup we focus on is the IIB duality frame with 03- and 07-planes with generalised 



fluxes defining a locally geometric background, i.e. F$, H, Q and P fluxes. In ref. 21 it was 
recently argued that, whenever one includes a number of fluxes equal to twice the number N 
of real fields in the theory, there is just enough freedom for casting the equations of motion 
into the form of a set of linear conditions, when the scalars at taken at the origin $o of 
moduli space. The space of solutions of the full problem will then have dimension N. The 
general solution will give the fluxes as a function of the N supersymmetry- (SUSY-) breaking 
parameters. 
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Table 1: Mapping between fluxes and couplings in the superpotential both in type IIB with 
03 and 07 and in type IIA with 06. The six internal directions of T 6 are split into " — " 
labelled by m = 1,3,5 and " | " labelled by a = 2, 4, 6. Note that the empty boxes in type IIA 
are related to the presence of dual fluxes which do not even admit any local description. Note 
that the orbifold involution forces i,j, k to be all different any time they appear as indices of 
fields of the same type (T or U). 

The SUSY-breaking parameters are N real constants A a and B a defined through 

D a W\* = A a + iB a . (6) 
In our non-isotropic locally geometric setup with N = 14 real fields, one would then need 



to consider 28 fluxes in order to apply the prescription of ref. 21 to the search for stable 
dS vacua. The set of 28 superpotential couplings that we have chosen here corresponds to 
the most general F% and H fluxes plus the first half of Q flux components in type IIB. From 
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table [TJ this set reads 



do, a l , a 2 , 03, Oo, o 1 , o 2 , 03, c , c x 



(7) 



The goal of this paper is to perform a scan of the 14-dimensional parameter space of 



solutions searching for fully stable dS extrema. In ref. 21 it was already argued that a 
random scan would not be efficient enough for finding special tachyon-free dS extrema, since 
the fraction of such points in the dS region is expected to be extremely tiny 

To overcome this difficulty, we have designed an evolutionary algorithm, which has the 
property of flowing in parameter space towards better-behaved solutions, i.e. positive cos- 
mological constant and positive eigenvalues of the mass matrix. This algorithm is presented 
in detail in appendix |B} 

Letting the algorithm run with the input parameters given in appendix [B] gave rise to 
3 fully stable dS vacua, which are shown in table [2j The full set of flux numbers for these 
solutions are specified in table [3] in appendix [A] . 
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1.85148 x 10 6 1.78128 x 10 6 
1.31064 x 10 6 1.27212 x 10 6 
113890 94187.9 
23907.3 11397.5 
5290.32 1478.37 
1353.35 607.799 
17.9045 2.85612 x 10~ 3 


262778 259704 
160140 153601 
24219.8 11296.3 
9273.53 7282.96 
4155.21 1745.39 
1306.57 343.391 
24.7202 11.4504 


34702.4 28731.5 

25215.5 18728.6 
5370.90 3609.62 
1572.22 1179.86 
723.060 518.923 
188.606 145.959 
5.16471 9.52228 x 10~ 4 



Table 2: The physical quantities characterising the 3 stable dS extrema found through the 
algorithm presented in appendix \B[ The first row shows the values of the cosmological con- 
stant, the second row the normalised energy j, and finally the third row shows the full mass 
spectra normalised to the cosmological constant. 

Finally we would like to carry out a complementary analysis w.r.t. the search that we 
have just presented. The aim is to map out the region of stable dS inside parameter space in 
the IIB case previously introduced and studied. By means of such an analysis one can study 
how stable dS vacua organise themselves in the parameter space spanned by {A a , B a }. 
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For each of the 3 stable dS extrema, we have plotted the level curves of 7 = ^j^l an d of 
77 = Min j Eigenvalues (ijn 2 ) 1 j in the (As, B$) plane. Subsequently we found the regions 
where 7 > 1 and 77 > 0, corresponding to stable dS critical points. These results are depicted 
in figure [T] for the three solutions. 

Our present understanding of the plots suggests that stable dS vacua organise themselves 
into thin sheets in parameter space. This conclusion seems to be perfectly in line with what 
was found in ref. 21 in the isotropic case where only N = Q real fields were retained in the 
model. 



Conclusions 

We have in this work presented three de Sitter vacua which are stable in all of the N = 14 
considered non-isotropic directions. These de Sitter vacua were obtained using the framework 
of the STU-models, more specifically Z2 x Z2 orbifold compactifications with generalised 
fluxes. We used a setup containing N = 14 real fields, with N = 14 supersymmetry- 
breaking parameters and 2N = 28 fluxes; F 3 , H and non-geometric flux Q. This result 
shows that the addition of non-geometric flux provides enough freedom to allow for stable 
de Sitter solutions. 

The details of the de Sitter solutions, the size of the cosmological constant and the values 
of the masses can be found in tabled We have also included table [3] for the reader interested 
in the details of the values of the fluxes. These solutions were obtained via an evolutionary 
algorithm designed to seek out stable de Sitter regions. We have included a description of 
our implementation of this algorithm in appendix [Bj 



As expected by ref. 21 , these solutions are part of stable de Sitter regions that are 
organised into sheets or small regions - an overlap between regions with positive cosmological 
constant and stability - as can be seen in 2-dimensional slices of the N = 14 dimensional 
parameter space. We have included plots of these regions in the vicinity of the found stable 
de Sitter solutions in figure [1} The structure and occurrence of these overlapping regions is 
very fascinating and is something we plan to investigate further in the future. 
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A Flux values and figures 





Sol. 1 


Sol. 2 


Sol. 3 


A s B s 
A Tl B Tl 
At 2 B Tl 
At 3 B T3 
Aui B Vl 
Au 2 B v . 2 
Au 3 B Us 


-0.314252 -0.677296 

-0.145819 0.201488 
0.324296 0.200242 

-0.510199 -0.830019 
0.437796 -0.223916 
-0.10521 -0.716205 

-0.110376 -0.952221 


0.266259 0.0513389 
0.0113161 1.19326 
-0.374583 -0.318198 
0.0745176 0.814969 
-0.410472 0.56115 
-0.074538 0.276484 
-0.597945 0.903104 


-0.586674 -0.186345 
0.684213 0.846619 
0.540869 1.10608 
-0.373167 0.724434 
-0.265422 0.187632 
-0.118099 0.651658 
-0.297359 0.75668 


a b 


0.249727 1.26278 


-0.0809612 -0.470688 


-0.126241 0.505546 


a« &« 


-0.385558 -1.00688 
3.24742 3.23091 
-3.42277 -1.65416 


0.763615 0.637699 
-0.129007 -1.32755 
-0.258342 -0.502742 


0.0845892 -0.342897 
-0.00376873 0.431055 
-0.687791 -1.81103 


af bf 


-0.678952 0.269654 
4.77196 -2.82032 
-4.74382 3.85504 


-1.69969 -0.585701 
0.141055 -0.0290134 
2.04863 0.623729 


-0.748501 -0.346495 
2.79635 -0.40546 
-2.13249 0.457822 




-0.626634 0.244097 


-0.0281293 -0.0877685 


-0.446200 -0.270327 




0.210031 
0.680146 
-0.154349 


0.224759 
-0.16114 
0.287961 


0.471208 
0.327864 
-0.586172 




0.934579 0.932678 0.134729 
1.3931 0.978188 0.351523 
-1.12313 -0.707561 -0.313208 


1.00141 -1.58335 -2.92823 
-0.335711 -0.460028 0.806899 
-0.628606 0.56901 1.78014 


-0.108522 0.144519 -1.70585 
0.178853 -0.083869 1.47008 
-0.489907 -0.220766 -0.305991 



Table 3: The 3 stable dS extrema found through the algorithm presented in appendix \B[ The 
first part of the table shows the values of the 14 SUSY-breaking parameters, whereas in the 
second part we give the explicit values of the 28 fluxes turned on at the corresponding critical 
points. 
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Figure 1: The parameter space of solutions around Sol. 1, 2 and 3 (top, center and bottom), 
projected on the (As, B$) plane. The origin in each picture is the found stable dS. Left: Level 
curves of the cosmological constant; the dS regions are the ones filled with lighter colours next 
to the Minkowski (blue dashed) lines. Middle: Level curves of the rj parameter; the tachyon- 
free regions are filled with darker colours. Right: The tiny region of overlap in parameter 
space corresponding with stable dS is zoomed in here. 



B The algorithm 



The search algorithm used is inspired by the methods of differential evolution. These are 
evolutionary algorithms that can walk in a solution landscape to find extremal points or 
special regions. Here we will describe the how the algorithm works and parts of the workings 
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of its implementation. We will start by presenting the steps of the algorithm, then proceed 
to explain them in closer detail. 

The algorithm proceeds according to the following repeated scheme. 

1. Generate a population (set of parents). 

2. Create a mutated copy of the population (set of children). 

3. Generate 7 and the mass eigenvalues for parents and children. 

4. For each parent and child pair, choose the better one. 

5. Check additional criteria and repeat until termination conditions. 

Step 1. The system is determined when the 14 SUSY-breaking parameters are specified. 
The population is a set where each element contain values for each of the 14 SUSY-breaking 
parameters. This set can either be generated randomly, or by a particular choice. We chose 
to leave one parameter to be fixed by demanding that we start at a Minkowski point. The 
size of the population is a weigh-off against the amount of iterations one wants to get to in 
a certain time. A large population gives a higher chance of a better initial position in the 
parameter space, but more iterations would more likely converge into a de Sitter region. 

Step 2. Performing a mutation is to take a step in some direction. This direction is 
random and so is the step size. This gives the algorithm the possibility to search freely 
without us imposing any explicit directions or step sizes. However, the randomness of the 
mutation should be chosen carefully - a mutation too random will make the algorithm a 
random-walk instead of being evolutionarily progressed. This is why the following code was 
chosen to mutate a member of the population 

mutate [memberO_ , sigmaO_] := Module [{member=memberO , sigma=sigmaO} , 
Do[ 

If [Not[RandomInteger[{l, 3}] == 1], 

member [[i]] = member [[i]] + RandomVariate [NormalDistribution [0 , sigma]] 

]; 

Length [member] }] ; 
member 

]; 

The step direction is random in the sense that there is a two in three chance that a number 
will be mutated. The size of the step is also random according to a normal distribution with 
zero mean. The standard deviation is given as input to this module, so that the algorithm 
can for example increase the effective step size if the solution gets stuck. There are hence 
two parameters we need to give to the algorithm that will affect its efficiency, i. e. the chance 
of getting mutated, and standard deviation for the mutation. 

Step 3. For each member we need the 7 and the mass matrix, to later be able to compare 
solutions. A module is needed to solve the set of equations and produce these parameters. 
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This is where most of the computational time is spent; given the SUSY-breaking parameters 
one needs to solve N linear equations and compute the eigenvalues oi a N x N matrix. 
This means that it is preferable to associate the 7 and each mass with every member of the 
population, to avoid calculating these more than once. 

Step 4. Some function must be invented to decide whether the parent or child is the 
better solution. The choice made here was to have the following three criteria 



where J\f(/i, a) denotes a random number chosen according to the normal distribution around 
the mean ji with standard deviation a. This applies to each of the 15 parameters that specify 
a stable dS extremum when all within range, i.e. 7 plus the eigenvalues of the mass matrix. 
All of the above conditions receive a number, its magnitude depending on how desirable the 
condition is. The total value parametrising the solution's desirability is given by the algebraic 
sum of the 15 separate contributions. This parameter will judge whether or not the parent 
should be replaced by its mutated child. If one of the above criteria is exclusively not satisfied, 
it receives the corresponding negative number instead. In choosing this evaluation criterion, 
some limit cases should be considered - what to do in the case of one parameter approaching 
range and one getting further away, etc.. Our choice of the normal distribution with mean 
A^ p , where the power p categorises the conditions by desirability is to make random choices 
in such conflicting cases. This is hence another parameter we need to give as input to the 
algorithm. 

Step 5. There are also additional conditions that could be checked for the code being 
able to progress smoothly for an indefinite time or until max iterations. Such a condition 
could be how long a certain solution have been iterated without finding a better child. If 
the solution has stayed in the same position for too long, one might want to increase the 
size of the mutation (the standard deviation), to jump out of this position where it seemed 
not able to find stable de Sitter. This gives two parameters to input into the algorithm: 
how long before a solution is categorized as stuck and how much the standard deviation 
should increase. One could also imagine that if the solution is very close of having 7 > 1 
and all non-negative masses, that one would decrease the size of the mutation to zoom in on 
a smaller search area, however this was not implemented by us. 

An additional check is to keep track of how long after a solution got stuck this should be 
removed and replaced. Since the search area is increasing after getting stuck, the chance of 
finding a better solution within reasonable time decreases. These solutions can be saved but 
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removed from the population, to be able to study them later if that is desired and to ensure 
the progression of the population. This gives another parameter as input to the algorithm. 
Similarly one should remove any stable de Sitter found, from the population and save it, this 
to make the population continue with the possibility to find more stable de Sitters. 

A termination condition should be given for the loop. The number of iterations, together 
with the populations size, is a comfortable way to determine the time for the code to run. 

As mentioned before, the algorithm has a number of parameters that needs to be specified. 
What follows are the values for these numbers which successfully produced stable de Sitter 
solutions. 

Population size 16 
Standard deviation for mutation a rn = 0.01 

Chance to mutate a number two in three 

Condition desirability Af(N p , a) (a = 1) 

Iterations 50 000 

Iterations without change until stuck 500 
Growth of mutation standard deviation per iteration after stuck .5er m 
Iterations without change until removed from population 2000 
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